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Abstract. We study the problem of computing the free space ^ of a simple legged robot called 
the spider robot. The body of this robot is a single point and the legs are attached to the body. The 
robot is subject to two constraints: each leg has a maximal extension R (accessibility constraint) 
and the body of the robot must lie above the convex hull of its feet (stability constraint). Moreover, 
the robot can only put its feet on some regions, called the foothold regions. The free space is the 
^\ ' set of positions of the body of the robot such that there exists a set of accessible footholds for which 

^>>^ ^ the robot is stable. We present an efficient algorithm that computes in 0(n^ log n) time using 

, 0{n^a{n)) space for n discrete point footholds where a{n) is an extremely slowly growing function 

(a(n) ^ 3 for any practical value of n). We also present an algorithm for computing when the 
foothold regions are pairwise disjoint polygons with n edges in total. This algorithm computes in 
^1^' 0(n'^ag(n) logn) time using 0{n^as{n)) space (ag(n) is also an extremely slowly growing function). 

O ' These results are close to optimal since Q{n'^) is a lower bound for the size of J^. 

o. 

■ 1. Introduction. Although legged robots have already been studied in robotics 
\ [ p^ n , only a very few papers consider the motion planning problem amidst obsta- 

O ■ cles ^ |7[ |l[ . In 1^, 1^ some heuristic approaches are described while, in Q efficient 

and provably correct geometric algorithms are described for a restricted type of legged 
robots, the so-called spider robots to be defined precisely below, and for finite sets of 
!>■ \ point footholds. 

A legged robot consists of a body with legs. Each leg has one end attached to the 
body and the other end (called the foot) that can lie on the ground (or move in space 

■ between two positions on the ground). Compared to the classic piano movers problem, 
&^ \ legged robots introduce new types of constraints. We assume that the environment 

' consists of regions in the plane, called foothold regions, where the robot can safely 

put its feet. A foothold is a point in a foothold region. The legged robot must satisfy 



o 
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c/3 ■ two different constraints: the accessibility and the stability constraints. A foothold 

^ ' is said to be accessible from a placement (position of the body of the robot) if it can 

be reached by a leg of the robot. A placement is called stable if there exist accessible 
k>( \ footholds and if the center of mass of the robot lies above the convex hull of these 

}_j ■ accessible footholds. The set of stable placements is clearly relevant for planning the 

I motion of a legged robot: we call this set the free space of the legged robot. Note that 

a legged robot has at least four legs, three legs ensure the stability of a placement and 
a fourth leg permits the motion of the robot. 



A first simple instance of a legged robot is the spider robot (see Figure 1.1). The 
spider robot was inspired by Ambler, developed at Carnegie Mellon University 
The body of the spider robot is a single point in the Euclidean plane and all its legs 
are attached to the body. The legs are retractable and their lengths may vary between 
and a constant R. We also assume that the center of mass of the robot is its body. 
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authorized footholds accessible footholds 




Fig. 1.1. The spider robot. 



It follows that a placement is stable if the body of the robot lies above the convex 
hull of the accessible footholds. 

The constraint that the body of the spider robot lies in the plane (instead of in 
3D) is not really restrictive. Indeed, consider a legged robot for which that constraint 
is relaxed. Then, if a placement {x,y,z) of such a legged robot is stable then, any 
placement {x,y,z'), ^ z' ^ z is also stable. Reciprocally, it can be shown that if 
{x,y) is in the interior of the free space of the spider robot, then there exists z > 
such that {x, y, z) is a stable placement of the corresponding legged robot. 

The problem of planning the motion of a spider robot has already been studied 
by Boissonnat et al. However, their method assumes that the set of footholds is 
a finite set of points and cannot be generalized to more complex environments. This 
paper proposes a new method for computing the free space of a spider robot in the 
presence of polygonal foothold regions. This method is based on a transformation 
between this problem and the problem of moving a half-disk amidst obstacles. Our 
method requires the computation of some parts of the free space of the half-disk. 
These computations are rather technical and complicated. Consequently, for the sake 
of clarity, we first present our algorithm for the simple case of discrete footholds, then 
we show how it can be generalized to the case of polygonal foothold regions. 

Once the free space of the spider robot has been computed, it can be used to find 
trajectories and sequences of legs assignments allowing the robot to move from one 
point to another. Indeed, once the free space is known, a trajectory of the body can 
be found in the free space. Then, a sequence of legs assignments can be computed as 
follows (sec [0 for details). Given an initial legs assignment, the body of the robot 
moves along its trajectory until it crosses the convex hull of its (three) feet that are 
on the ground or one leg reaches its maximal extension. Then, a suitable foothold is 
found for the fourth leg and one leg leaves its foothold. 

The paper is organized as follows: some notations and results of |Q| are recalled in 
the next section. Section ^ shows the transformation between the spider robot problem 
and the half-disk problem. We present in Section ^ our algorithm for computing the 
free space of a spider robot for a discrete set of footholds. Section [H] shows how to 
extend the algorithm to polygonal foothold regions. 

2. Notations and previous results. In Sections 2, 3 and 4, S denotes a dis- 
crete set of distinct footholds {si, . . . ,Sn} in the Euclidean plane {S will denote in 
Section 5 a set of disjoint polygonal regions). Point G denotes the body of the robot 
(in the same plane) and [0, R] is the length range of each leg. The free space !F is the 
set of all stable placements of G. A placement is said to be at the limit of stability if 
it lies on the boundary of the convex hull of its accessible footholds. Notice that is 
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a closed set and contains the placements at the limit of stability. 

Let Ci denote the circle of radius R centered at Si . A is the arrangement of the 
circles Ci for 1 ^ i ^ n, i.e., the subdivision of the plane induced by the circles. 
This arrangement plays an important role in our problem and we will express the 
complexity results in term of \A\, the size of A. In the worst-case, \A\ = 8(n^) but if 
k denotes the maximum number of disks that can cover a point of the plane, among 
the disks of radius R centered at the Si, it can be shown that |^| = 0{kn) p^ . 
Clearly k is not larger than n and in case of sparse footholds, |^| may be linearly 
related to the number of footholds. 

For any set £, let d{£) denote its boundary, CH(f ) its convex hull, int{£) its 
relative interioi]^ clos{£) its closure, and compl{£) its complementary set. Let 
denote the set of angles ]R/2ttZ. We denote by x = y[p] the equality of x and y 
modulo p. We say in the sequel that two objects properly intersect if and only if their 
relative interiors intersect. 

The algorithm described in [0 is based on the following observation: for G in a 
cell F of A, the set of footholds that can be reached by the robot is fixed; the portion 
of F that belongs to !F is exactly the intersection of F with the convex hull of the 
footholds that can be reached from F. Therefore, the edges of d{T) are either circular 
arcs belonging to A or portions of line segments joining two footholds. Moreover, 



a vertex of d{J-') incident to two straight edges is a foothold (see Figure 2.1). The 
complexity of T has been proved to be \T\ = 6(|-4|) 




Fig. 2.1. An example of the free space of a spider robot. 

The algorithm presented in jl] computes the free space J- in 0(|yl| logn) time. It 
uses sophisticated data structures allowing the off-line maintenance of convex hulls. 

The algorithm described in this paper has the same time complexity, uses simple 
data structures and can be extended to the case where the set S of footholds is a set 



^The relative interior of a set £ in a space E is the interior of £ in the space £ for the topology 
induced by E. For example, the relative interior of a closed line segment in is the line segment 
without its endpoints, though its interior in IR^ is empty. 
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of polygonal regions and not simply a set of points. For simplicity, we consider first 
the case of point footholds and postpone the discussion on polygonal foothold regions 
to Section ^. 

General position assumption. To simplify the presentation of this paper, we 
make the following general position assumptions. All these hypotheses can be removed 
by a careful analysis. Recall that we consider here that the set of footholds is discrete. 

No two footholds lie at distance exactly R or 2R. Among the circles Ci , . . . , C„ 
and the line segments joining two footholds, the intersection between three circles or, 
two circles and a line segment or, one circle and two line segments, is empty. 

3. From spider robots to half-disk robots. In this section, we establish the 
connection between the free space of the spider robot and the free space of a half-disk 
robot moving by translation and rotation amidst n point obstacles. 

Theorem 3.1. The spider robot does not admit a stable placement at point P if 
and only if there exists a half-disk (of radius R) centered at P that does not contain 



any foothold of S ( see Figure 3.1) 



Proof. Let TZ be the set of all the footholds that are reachable from placement 
P. By definition, P is not stable if and only if the convex hull of TZ does not contain 



P (see Figure 3.1). That is equivalent to say that there exists an open half- plane 
through P containing TZ, or that there exists a closed half-disk of radius R centered 
at P which does not contain any foothold. □ 



Fig. 3.1. A placement which is not stable. 




Fig. 3.2. HD{P,e). 



Definition 3.2. Let HD{P,9) be the half-disk of radius R centered at P (see 
Figure \3.^) defined by: 



{x-xpY + {y~ypf ^R^ 

{x — Xp) sa\9 — {y — yp) cos9 ^ 
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Definition 3.3. Vs^ £ S {1 ^ i ^ n) let: 

n 

i=l 

C i = Ci X S . 



Tii will be called the helicoidal volume centered at s,; (see Figure S.o). 



e 




Fig. 3.3. Helicoidal volume "Hi, 

Notice the typographical distinction between the circle Ci defined in and the 
torus Ci defined in x 5^. For convenience, we will often identify and the interval 
[0, 27r] of M. This allows us to draw objects of IR^ x 5*^ in IR^ and to speak of the 
6'-axis. Ilg^ denotes the "plane" {iP,0) £ x \ 9 ^ eo}- 

Definition 3.4. The free space C of a half-disk robot moving by translation and 
rotation amidst the set of obstacles S is the set of (P, 6) G x such that the 
half-disk HD{P, + tt) does not intersect S. 
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Proposition 3.5. L — complCH). 

Proof. ¥6* e S'\ the set £ n He is the free space of the half-disk HD{P, 6 + tt) 
moving by translation only, amidst the obstacle si, . . . , s„. Since the set of points P 
such that HD{P, 6 + n) contains a is HD{si, 0), C Olle is the complementary set 
of the union of the HD{si, 9). Thus, C is the complementary set of the union of the 
Hi, that is H. □ 

Let pi/g denote the mapping (called "orthogonal projection"): x — > 
]R^,{P,e) ^ P. 

Theorem 3.6. T = compl{p^/g{compl{H))) 

Proof. By definition of C, p//g{jC) is the set of points P € such that there 
exists an angle 6 G such that the half-disk HD{P, 0) does not intersect S. By 
Theorem 3.1, it is equivalent to say that there exists 9 E such that HD{P,9) 
does not intersect S, or that P is not a stable placement of the spider robot. Thus, 
p//g{£) is the set of points P where the robot does not admit a stable placement, i.e., 
= compl{p//g{£)). The result then follows from Proposition 3^. □ 

Remark 3.7. compl{p//g{compl{'H))) x is the largest "cylinder" included 
in H, whose axis is parallel to the 0-axis (in grey in Figure iA). The basis of this 
cylinder is J-. 




compl{p ^^g{compl {£))) 
Fig. 3.4. compl(pifg{compl{£))). 



Remark 3.8. The results of this section do not depend on the fact that the 
footholds are discrete points. For more general foothold regions, we simply need to 
replace the helicoidal volumes by their analog. This will be done in Section ^. 

4. Computation of T. In this section, we propose an algorithm for computing 



J- based on Theorem 3.6 



A first attempt to use Theorem 3^ may consist in computing C = compliTi) and 
projecting it onto the horizontal plane. The motion planning of a convex polygonal 
robot in a polygonal environment has been extensively studied (see for example jlo| . 
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Such algorithms can be generahzed to plan the motion of a half-disk. It should 
lead to an algorithm of complexity 0(nXs{n)logn), where Xs{n) is an almost linear 
function of n. The projection can be done using classical techniques, such as projecting 
all the faces of C and computing their union. Since the complexity of the 3D object 
C is not directly related to the complexity of its projection, this approach does not 
provide a combinatorial bound on JF. However, assuming \J-'\ = 0{Xs{\A\)) (which 
will be proved in this paper) the time complexity of the algorithm of Kedem et al. is 
0{nXs{n) logn + A^(l^l)log^n). 

In this paper, we present a direct computation of J-'. This approach provides an 
upper bound on the size of J-', namely \ J- \ = 0{Xs{\A )). It also provides an algorithm 
for computing T in 0(As(|^|) logn) time. As in and contrary to ||ll[, the algo- 
rithm proposed here is sensitive to |^| which is usually less than quadratic. Another 
advantage of our direct computation is to avoid the explicit construction of the 3D 
object £ which is useless for our application. Our algorithm manipulates only two- 
dimensional arrangements or lower envelopes and we provide a detailed description of 
the curves involved in the construction. 

Let us now detail the computation of in the case of point footholds. We know 
that each arc of the boundary d{J^) of ^ is either a straight line segment belonging 
to a line joining two footholds or an arc of a circle Ci (see Section ||). The circular 
arcs d{J-') D Ci are computed first (Sections |4.1|, [4.21 and 



the line segments in a second step (Sections 4.4 



and 4.5) 



4.3) and linked together with 



4.1. Computation of d{T) A. In the sequel, the contribution of an object 
X to another object Y is X DY. We compute the contribution of each circle Cj^ , 
io = to d{J^) in turn. Recall that Cig denote the torus Cig x S^. The 

contribution of each circle Cig to d{T) will be obtained by computing the intersection 
of all the Tii, i — 1, . . . ,n, with the torus Cig. Let Zi, i — 1, . . . denote these 
intersections: 

2i = ^M-i n Cig . 

We first show how to compute the contribution of Cig to d{J-') in term of the Zi, 



and leave the studies of the shape and properties of Zi to Section 4.2. Figures tllJ 



and 4.2 show some (hatched) Zi C Cig {i ^ io) where Cig is parameterized by (m, 0) 
(u and 9 parameterize Cig and respectively); the dark grey region shows Zig. 
Proposition 4.1. The contribution of Cig to d{T) is: 

Cig n d{T) = compl{p//Q{compl(UiZi))) \ int{compl{p//g{compl{Ui^igZi)))) . 



Proof. Since is a closed set, Cig n 9(JF) — [Cig n JF] \ [Cig n int{T)]. According 
to Theorem ^.6| , T — compl(p//Q{compl{Ti.))). One can easily prove that for any 
set £ G X , int{compl{£)) = compl{clos(£)), clos{compl{£)) = compl{int{£)), 
and clos{pi/g{£)) = p / j g{clos{£)) . It then follows from the expression of J- that 
int{^F) = compl{p/^g (comprint (Ti.)))). 

Recall that for any sets X,Y € IR^ x S^, compl{X n F) = compl{X) U compliY), 
p//g{X U y) = p//g{X) U P//g{Y), and compl{X U F) = compl{X) n compl{Y). That 
implies 

compl{p//g{compl{X n Y))) — compl{p //g{compl{X))) n co7npl{p //g{compl{Y))). 
We now consider that equation with X equal to Ti or int{'H), and Y equal to the 
torus Cig. Since compl{p//g{compl(Cig))) is the circle Cig we get: 
compl{p^/g{compl{H n Cig))) =J^r\Cig and 
compl{p^/g{compl{int{H) H Cig))) — int{T) n Cig. 
Since U = iJ'"^l'H^ and Z, = n C^g by definition, H n C^g = Uf^^Z, and 
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int{T-L) n Cig — U"^^{int(T-li) n Ci„). By the general position assumption, no two 
footholds lie at distance 2R, thus for i 7^ iq, int{'Hi) OCig = mt(Zi)0. As int{'Hig) D 
Cig = 0, we get int{Tl) n Cig — Ui^i„int{Zi). The study of the shape of Zi will 
yield (see Lemma [4.8[ ) that Ui^igint{Zi) = int{\Ji:^igZi). Therefore, int{!F) H Ci^ = 
compl{p//g{compl{int{Ut^iaZi)))) = int{compl{p//B{compl{\Ji^i^Zi)))) and J^nCio = 
compl{p//g{compl{LliZi))). Using Ci„ n d{J^) — [Ci„ n JF] \ [Ci^ n int{J^)], we get the 
result. □ 

Thus, the contribution of dg to d{J-) comes from the computation of UiZi and 
Geometrically, compl{jp//g[compl[yJiZi))) is the vertical projection (along the 6- 



axis) of the largest vertical strip Ej^ included in UiZi (see Figure 4.2). Similarly, 
compl{p//g{compl(Ui^igZi))) is the projection of the largest vertical strip H'^^ included 
in Ui^igZi. Thus, d{J-') D Ci^ is the vertical projection onto Cig of the vertical strip 
S,„\znt(S^J. 

In order to compute J- efficiently, we need to compute the union of the regions 
Zi efficiently. More precisely, we will show that the union of the regions Zi can be 
computed in 0{kig\ogkig) time where kig is the number of helicoidal volumes TLi 
intersecting Cig . 

This is possible because the Zi have special shapes that allow us to reduce the 
computation of their union to the computation of a small number of lower envelopes 
of curves drawn on Cig , with the property that two of them intersect at most once. 



The geometric properties of the Zi are discussed in Section and, in Section |4.3| , 
we present and analyze the algorithm for constructing d{!F) n Cig. 

4.2. Properties of the Zi. We study here the regions Zi = TLi C^Cig. Recall 
that we parameterize Cig = Cig x by (w, 9) where u and 9 parameterize dg and 
resp ectiv ely (u = corresponds to the point of dg with maximum x-coordinate) . 



Figures and [4.2| show examples of such regions Zi. For convenience, we will use 
the vocabulary of the plane when describing objects on the torus Cig. For instance, 
the curve drawn on the torus Cig with equation a9 + bu + c = will be called a 
line. The line u = uq will be called vertical and oriented according to increasing 9. 
Lower and upper will refer to this orientation. The discussion below considers only 
non empty regions Zi (such that ||sipSi|| < 2R). 

We introduce first some notations. Let HCi{9) be the half-circle of the boundary 
of HD{s^,9), i.e., HC^{9) ^C^C^ HD{.%,9). Let ri{9) be the spoke of that makes 
an angle 9 with the x-axis, i.e., ri{9) — {si + Xug \ A G [0, i?]} where ug is the unit 
vector whose polar angle is 9. The boundary of Tii is composed of the three following 
patches: 

% = {{HCm,9)^R'xS^} 
T^t = {{n{9),9) eM^x S^} 
m = {{n{9 + 7r),9) eM^x S^} 

Let p~ and pf denote the curves TZ^ n dg and TZf n dg , respectively. Since TZ~ 
and TZf are translated copies of one another, i.e., TZ~ = TZf ± (0, 0, tt), we have: 

Lemma 4.2. p~ and p^ are translated copies of one another, i.e., 
pt = {{u, 0) e 51 X 51 1 (u, 9-tt)^ p-} = {{u, 0) e 51 X 51 1 [u, 9 + Tr)e Pi}. 



^Recall that int denotes the relative interior, thus int{'Hi) is the interior of Hi in IB? X but 
int{Zi) denotes the interior of Zi in dg. 
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Lemma 4.3. The curves pf are monotone in u. 

Proof. Assume for a contradiction that a curve is not monotone in u. Then, 
there exists u and 9 ^ 6' in such that {u,9) and {u,9') parameterize points of pf . 
By the definition of TZf , it then follows that the point U G Ci^ parameterized by u 
belongs to the two spokes ri{9) (or ri{9+TT)) and ri{9') (or ri{9'+n)). The intersection 
between any two of these spokes is exactly Si. Thus, U = Si, which contradicts (since 
U S Ci„ ) the general position assumption saying that the distance between Si and Sig 
is not R. □ 

Lemma 4.4. The region Zt^ is the subset of C j^ p aram eterized by {{u,9) G 5^ x 
\ 9 ^ u ^ 9 + tt} (shown in grey in Figures 4^.1 and \4.^ ). 

Proof. For any 9 G S^, the intersection between Hig and the "horizontal plane" 
He is the half-disk HD{sig, 9). Similarly, the intersection between Ci„ and that plane 
is Cig. Thus, the intersection between Zi^ and Wg is HCig{9), which is parameterized 
on Cig hy {u G \ 9 ^ u ^ 9 + Tr}. That intersection is actually on the plane Ilg and 
is therefore parameterized on Ci^ hy {{u,9) E x \ 9 ^ u ^ 9 + n}. □ 

Proposition 4.5. Zi is a connected region bounded from below by p~ and from 
above by pf , i.e., Zi = {{u,9) € 5^x5^ | 3x S [0,7r], {u,9-x) E p'[ , {u,9-x+Tr) E pf} 
(see Figures 4jJ_, U-^ - 

Proof. By cutting Ci^ and Tii by the "horizontal plane" Ilg, we get that a point 
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parameterized by (w, 9) on Cig belongs to Tii if and only if the point U parameterized 
by u on Cig belongs to HD{si,9). Since HD{si,0) can be seen as the union of the 
spokes {ri{6 + 7) | 76 [0, tt]}, {u,6) e Zi if and only if there exists 7 G [0, tt] such 
that U e ri{9 + 7), or equivalently, U G ri{9 — x + tt) with x = tt — 7 g [0, tt]. Since 
7^7 = {(ri(6' - a; + tt), 6' - x) I 6* - x G S^}, it follows from C/ e ri(6' - x + tt) that the 
point of Cig parameterized by {u,9 — x) belongs to TZ~ and thus to = TZ~ n Cig. 
From Lemma we get that the point parameterized by {u,9 — x + tt) belongs to 
pf . Therefore, Zi is a connected region bounded from below by p^ and from above 
by p+. □ 

We want to compute the union of the Zi by computing the "lower envelope"]^ of the 
lower edges and the "upper envelope" of the upper edges pf . It is unfortunately 
impossible to do so because some upper edges pf may possibly be "below" or intersect 
some lower edges pj. However, we can subdivide the regions Zi into blocks Zi , k G IC, 
and separate these blocks into two sets fii and such that the union of the Zf in 
fli (resp. 512) is the region bounded from above by the upper envelope of the upper 
edges of the Zf G fli and bounded from below by the lower envelope of the lower 
edges of the G fli (resp. fl2)- Such property can be realized by showing that all 
the upper edges of the Zj^ G fli belong to the strip {{u,9) G x [u+ ^,u+ ^]} and 
all the lower edges of the Zf G fii belong to the strip {{u, 9) G x [u — ^,u + ^]} (a 
similar property is shown for 0.2)- Note that the upper and lower envelopes are then 
defined since they are considered in x M. 




Fig. 4.3. For the definition of r[{9) and r"(9). 



We subdivide Zi into blocks Zi when R < \\sigSi\\ < \/2R. That subdivision 
is performed such that the upper and lower edges of the are 6'-monotone. Recall 
that the upper edge pf of Zi is the intersection of TZf = {{ri{9),9) \ 9 G S^} and 
dg. The spoke ri{9) intersects Cig twice (for some 9) when R < \\sigSi\\ < V^R, 
which implies that pf is not 6'-monotone. We cut the spoke ri{9) into two pieces 
such that each piece intersects Cig at most once. Let T be the intersection point 
between Cig and on one of the two lines passing through Si and tangent to Cig (see 
Figure [4.3[ ). Let T(9) be the point on ri{9) at distance |jsiT|| from Si. Cutting 
ri{9) at T{9) defines two sub-spokes r^{9) and r'/{9) that intersect Cig in at most 
one point each; without loss of generality, let r[{9) denote the sub-spoke joining Si to 



^Note that the lower and upper envelopes of curves in X are not actually defined. 
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T{9). The set oi 6 E for which r'^{0) intersects Cig is clearly connected but the 
set oi 9 E for which r"(6') intersects Cig consists of two connected components. 
We denote by the intersection {{r[{e),e) \ 9 € S^} r\C^„ and by p-+ and pf+ 
the two connected components of the intersection {{r'l{9),9)\ 9 G S*^} n C^^ (see 
Figure 1^3). Since r[{9) and r'l{9) intersect at most once for any 9 & S^, the 
curves pf^ , pf^ and p^^ are 6'-monotone. The lower edges k = 1,2,3 can be 
defined similarly or in a simpler way as the translated copies of p'l'^ , k — 1,2, 3, i.e., 
pfe" = {{u, 9)€S^ xS^\ {u, 9 + n) e p+'^j. We denote by , fc = 1, 2, 3, the subset 
of Zi bounded from above by p^'^ and from below by p*^^ . 

We can now prove the following proposition that will allow us to compute the 
union of the Zi by computing the upper and lower envelopes of their upper and lower 
edges. 

Proposition 4.6. //O ^ llsio'^ill < line 9 ^u - ^ properly intersects Zi, 

and the lines 9 — uzt ^ properly intersect neither p'f nor p~ . 

If R < \\sigSi\\ < y/2R, the line 9 = u + ^ properly intersects Zf , and the line 
9 — u — ^ properly intersects Zj and Zf . Furthermore, the lines 9 = u± ^ properly 
intersect none of the edges p]^ , p\^ , pf^, p^^ , p^^ and p^^. 

If \/2 R ^ llsioSill < 2R, the line 9 ^ u + ^ properly intersects Zi, and the lines 
9 = u± ^ properly intersect neither pf nor p^ ■ 



Proof. Let {up,9p) parameterize a point of a curve pi 



Let P denote the point 
One can 



of Ci„ with parameter up and 7 = Z{Psi„, Psi) [27r] (see Figure |4.4D 
easily show that 7 = 9p — up[t:]. We prove that 77^ ^ [tt], except possibly when 
{up,9p) is an endpoint of p^ (or pf when R < \\sigSi\\ < \/2R), which implies, since 
7 = 0p — Mp[7r], that the lines 9 = uzL ^ intersect neither pf nor p~ (resp. 
Pi^), except possibly at their endpoints. 



p^'+ nor 





Fig. 4.4. For the proof of Proposition |^. fij 



Case 1: ^ |lsio'*j|l < ^- Since Si belongs to the disk of radius R centered 



at Sio, 7 e (- 



2 ' 2 . 



for any P e d, 



properly intersect neither nor p^ 



(see Figure AA). Thus, the lines 9 
Finally, the point of Cig (^2,^*2 — | 



^±1 
where 

and also to the relative interior of 



6*2 = Z{x, Si^i) [2tt], belongs to the hne 9 ~ u^ 

Zi since it belongs to the interior of Ti^ (see Figure 4.5 a). Therefore, the line 
properly intersects Zi. 

Case 2: R < \\sigSi\\ < \/2R. Let {up^,9p^) parameterize the point connecting 
p\^ and pf^ , and {up.^,9p2) parameterize the point connecting pf^ and p^^. Let Pi 
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(c): R < II s.qS, II < ^/2R, (01, di - f ) 6 Z} 



(d): R < ||s,qS,|| < s/2R, (63,63 - f ) G 2f 



Fig. 4.5. For the proof of Proposition 4.6: section of "Hi and dg by the "planes" Ug^ 



Ilgg _ ^ respectively. 



and P2 denote the points of Ci„ parameterized by u and up^ respectively. According 
to the construction of p^"*", pf^ and p'^'^ , the tangent hnes to Cig at Pi and P2 pass 
through Si. At most two tangent hnes to Ci^ pass through s^, thus Pi and P2 are 
the only points of Ci„ where 7 = f [tt]. Since is u-monotone by Lemma 4.3 , 
[up^,9p^) and {up.^,6p.^) are the only points of pf where 7 = f [tt]. Therefore, the 
lines 9 = u±^ do not properly intersect , fc = 1,2,3. Similarly, the lines 9 — uzL^ 
do not properly intersect p^~, fc = 1, 2, 3. 

Let 01 and 9^ be the parameters on of the intersection points between dg 
and Ci (see Figures 4.5c and d); to differentiate 9i from ^3, assume without loss of 
generality that, for any e > small enough, the points of dg parameterized by 6'i +e 
and 03 — £ are in the disk of radius R centered at Si. Then, the points (0i, 6*1 — ^) and 
'^3, ^3 — f ) of Cig belong to Zf and Zf (or to Zf and Zf) respectively (see Figures 4.5 



and d). However, these points do not belong to the relative interior of Zf and Zf 
(because they lie on the border of HD{si, 9i — ^) and HD{si,9^ — f )). Nevertheless, 
there clearly exists e > small enough such that the point parameterized by 6*1 + e 
(resp. 6*3 — e) on Cig belongs to the interior of the half-disk HD{si, 9i — ^ + e) (resp. 
HD{si, 6*3 - f - e)). Thus, the points (0i + e, 6*1 + e - f ) and (6*3 - e,\ - e - f ) 
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of Cig belong to the relative interior of and Zf respectively. Therefore, the line 
9 ^ u — ^ properly intersects Z^ and Zf. 

On the other hand, (6*2, 6*2 + f) (where 02 = Z{x,Si^i) [2tt]) belongs to relative 
interior of Zf because the point of Cig parameterized by 02 belongs to the relative 
interior of the sub-spoke r[{92 + tt) (see Figure p^b) which belongs to interior of 
HD{si,62 + f ). Therefore, the line 9 = u + ^ properly intersects Zf. 

Case 3: ^/2R ^ || s^oSiH < 2R. Since ri{9) intersects Cig at most once, 76 [f , 
(see Figure 4^). Moreover, 7 = f [tt] only when ||si(,Si|| — V2R, but then, P is at 
distance R from Si which implies that {up,9p) is an endpoint of pi. Thus, the lines 
9 = intersect neither nor p~ , except possibly at their endpoints. Finally, the 

point (6*2, 6*2 + f ) of Cig (where 6*2 = ■^{x, sj^i) [2 7r]) belongs to the hne 9 — u + ^ 
and also to the relative interior of Zi (see Figures 4.5b and 1.1). Therefore, the line 
9 = u + ^ properly intersects Zi . □ 

By Proposition 4.6, we can compute the union Ui^igZi by separating the Zi, 
Z^ into two sets Qi and Q2 (where Zi, Zf belongs to fii if and only if 
belongs to the strip {(-u, 0) S S*^ x [m + -I, m + ^]} and p~ , p^~ belongs to the strip 
{{u,9) e S*^ X [u — ^, u+ ■!]}) and computing the union of the Zi, Z\ in fii (resp. ^2) 
by computing the upper envelope of their upper edges and the lower envelope of their 
lower edges. In order to compute efficiently these upper and lower envelopes, we show 
that the curves p'l , p^ , p^'^ and p\~ intersect each other at most once. However, we 
need for that purpose to split the regions Zi when < HsigSiH < R into two blocks 
Z} and Z'^ separated by the vertical line u — 92 — Z{x, Sigs]); it also remains to split 
the ^-interval (or the it-interval) over which pi is defined into two intervals of equal 
length over which p\'^ and pf^ are defined (see Figure 4.2). Note that Proposition |4.6| 
still holds if we replace (when < ||sioSi|| < R) Zi hyZf and pf by p*^^, fc = 1,2. 

For consistency, we split Zig into two blocks Zl^ and Zf^ separated by a vertical 
line (chosen arbitrarily, say u = n). Also for consistency, the curves pf^ when \/2 R ^ 
||si(,Si|| < 2R are occasionally denoted in the sequel pl^ . 

Lemma 4.7. Let p[ and p'j be some connected portions of p^ and p^ respectively 
fi 7^ j)- U Pi or p'j is monotone in 9 and defined over a 9-interval smaller than n, 
then p'j^ and p^ intersect at most once. 

Proof. Let {ui, 9j) be a point of intersection between p,^ and p'j and / be the point 
of the circle Cig with parameter uj. Since p^ is a portion of the intersection between 
Cig and TZf , / is a point of intersection between Cig and the diameter of HD{si, 9i). 
Therefore, the line passing through Si and / has slope 9j. 

By applying the same argument to pj, we obtain that Si and Sj belong to the 
same straight line of slope 9i. Therefore, if p^ and pj intersect twice, at {ui,9i) and 
{uj,9j), then 9i = 9j[k\. It follows, if p[ or p^ is defined over a 0-interval smaller 
than TT, that 9i ~ 9j[2'!i\. Furthermore, if pj or p^ is monotone in 9, then (u/, 9i) and 
{uj,9i) are equal. □ 

Lemma 4.8. Vi, j, int{Zi) U int{Zj) = mt(Zj U Zj). 

Proof. We assume that i ^ j because otherwise the result is trivial. One can 
easily show that int[Zi) U int{Zj) ^ int{Zi U Zj) only if the boundaries of Zi and Zj 
partially coincide, i.e., the dimension of d{Zi) D d{Zj) is 1. 

By Proposition |4.5| , d{Zi) consists of the edges pf and p~ and of two vertical line 
segments joining the endpoints of pf and p~ when these endpoints exist (which is the 
case when i ^ iq). Moreover, these vertical line segments are clearly supported by 
the vertical lines u — 9i and u — 9^ where 9i and ^3 parameterize on Cig the points 
of intersection between Cig and Ci (see Figure 4.1). 
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By Lemma 4.7, the edges pf and pj do not partially coincide. By the general 
position assumption, no three distinct circles C^q, Ci and Cj have a common inter- 
section point. Thus, for any i ^ j, D Ci and Ci^ n Cj are disjoint. Therefore, the 
vertical hues d{Zi) \ {p'l , p~} and d{Zj) \ {p'^ , pJ} do not partially coincide. Finally, 
since pf is nowhere partially supported by a vertical line by Lemma 4.2, pf and the 
vertical lines d{Zj) \ {p^ , pJ} do not partially coincide. □ 

Proposition 4.9. Any two curves among the curves p^^ intersect at most once 
(where k e {1,2} i/0 ||sioS»ll < R, k e {1,2,3} if R < \\s^„s^\\ < V^R, and k = 1 
ifV2R^\\s,,s^\\ <2R). 

Proof. By Lemma |477| , it is sufficient to prove that all the curves p*'^, i 7^ «o, are 
monotone in 9 and defined over 0-intervals smaller than tt. Indeed, the curves p}^ , 
Pio ^ Pit ^'^^ Pio clearly do not pairwise intersect more than once, by Lemma 4.4. 

If < ||si(,Si|| < i?, any spoke of d intersects at most once. Hence, pf is 
monotone in 9. pf is defined over a 0-interval greater than tt but smaller than 2it. 
Since we have split that interval in two equal parts, p\^ and p1^ are defined over a 
^-interval smaller than tt (see Zi in Figure 4.2). 

li R < \\si„Si\\ < v^i?, the ^-interval where ri[9) (or ri{9 + vr)) intersects Ci„ is 
smaller than tt, which implies that pi is defined over a 0-interval smaller than tt. The 
curves p^^ ^ k = 1,2,3, are defined as the connected components of {{r[{9),9) \ 9 G 
S'^}C^Cia and {{r'l{9),9) \ 9 e S^}C^C^„. Since the sub-spokes r[{9) and rf (6*) intersect 

at most once for any 9 lE S^, the curves p*^^, fc = 1,2,3, are 0-monotone. 

If V^R ^ ||si(,Si|| < 2i?, ri{9) (and also ri{9 + tt)) intersects Ci„ in at most one 
point, which proves that pi is monotone in 9. Furthermore, the (^-interval where pi is 
defined is smaller than tt because the 6'-interval where ri{9) (or ri{9 + tt)) intersects 
Ci„ is smaller than tt. □ 

4.3. Construction of d{f^) D C,„. We first show how to compute UiZi. Let fii 
and be the following sets of Z,f : 

ill = {Z, I V2R < ||s,„s,|| < 2R} U {Zf I R < \\s,„s,\\ < V2R}, 
^2 = {Zl,Zf I ^ \\s.,s,\\ <R}U {Zl,Zf I R < \M\ < V2R}. 

By Proposition |4.6| , the line 9 = properly intersects all the G fli but the 

lines 9 = u± ^ properly intersect none of their upper and lower edges p^"*" and . 
Thus, the regions Z^ e f^i can be seen as regions oi {{u,9) e 5^ x [u ~ + ^]} 
such that all their upper edges p'^'^ he in {{u,9) £ x [u -f- -I, u -t- ^]} and all their 
lower edges p^~ lie in {{u,9) € x [u — f ,u + ^]}. Therefore, the union of the 
Z^ E ill is the region of {(it, 0) e 5^ x [u — + ^]} bounded from above by the 
upper envelope of their p'l^ and bounded from below by the lower envelope of their 
p'^~. Similarly, the union of the Z^ g ^2 is the region of {(u, 6*) e 5^ x [u - ^, u-f- f ]} 
bounded from above by the upper envelope of the p'l^ and bounded from below by 
the lower envelope of the p^~ . 

The union of Oi and Q2, which is UiZi, can be achieved by computing, on one 
hand, the intersection between the upper edge chain oiU^i'^ci-^Z^ with the lower edge 

chain of U^^g^^^Zf (which both belong to {{u,9) e x \ 9 e [u + + ^]}), 
and on the other hand, the intersection between the upper edge chain of ^zf'en2-^i 
with the lower edge chain of ^z''eni^i (which both belong to {(m, 9) & x \ 9 £ 
[u — j,u + ^]}). These intersections can simply be performed by following the two 
edge chains for u from to 27r, since they are monotone in u by Lemma [4.3| . 

Let us analyze the complexity of the above construction. The ki^ helicoidal vol- 
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umes TLi that intersect C^q can be found in 0{ki^) amortized time once the Delaunay 
triangulation of the footholds has been computed, which can be done in O(nlogn) 
time |, |l|l. By Proposition U, the upper and lower envelopes can be computed 
in 0{kig logfcig) time using 0{kiga{kig)) space where a is the pseudo inverse of the 
Ackerman's function ||]. Also by Proposition 49, the union of ili and can be 
done in linear time in the size of the edge chains, that is 0{kiga{kig)) time. Thus, 
we can compute UiZi in 0{ki„ logfei^) time using 0{kiga{kig)) space after 0(nlog?T.) 
preprocessing time. We can compute Ui^i„Zi similarly by removing Z^^ and Zf^ from 



The contribution of Ci^ to d{J-) is, according to Proposition 4.1, Ci^ Qd{J-) 



C(jmpl{p / 1 g{ccympl{yJiZi))) \ int{ccympl(j) / j g{c(jmpl{\Ji^igZi)))) . By Remark 3.7 



compl(jj//g{compl(UiZi))) and compl{p//g{compl(Ui^igZi))) are the projections onto 
Cig of the largest vertical strips and included in UiZi and Ui^i^Zi, respectively 
(see Figure 4.2). These projections are easily computed because the edges of UiZi 
and Ui:^igZi are monotone with respect to u (Lemma |4.3| ). These projections, and 
therefore the computation of Ci^ D d{J^), can thus be done in linear time and space 
in the size of UiZi and Ui^i^Zi, that is 0{kiga{kig)). 

Moreover, we label an arc of d{!F) either by i if the arc belongs to the circle Ci or 
by («, j) if the arc belongs to the straight hue segment [si, Sj]. The labels of the edges 
of d{J^) incident to Cig can be found as follows, without increasing the complexity. 
An arc of d{J^) (1 Ci^ corresponds to a vertical strip T,ig \ S . An endpoint P of such 
an arc is the projection of a vertical edge, or the projection of a point of intersection 
between two curved edges. In the first case, P is the intersection of Cig with some Ci 
and in the second case, P is the intersection of Cig with some line segment [s^, Sj]. By 
the general position assumption, among the circles Ci , . . . , C„ and the line segments 
joining two footholds, the intersection between three circles or, two circles and a line 
segment or, one circle and two line segments, is empty. Thus, P is the intersection 
between Cig and either a unique Ci or a unique line segment [si, Sj]. Therefore, the 
edge of d{T) incident to dg at P is either a circular arc supported by d or a line 
segment supported by [si, Sj]. Hence, the labels of the edges of d{!F) incident to Cig 
can be found at no extra-cost during the construction. 

Since A is the arrangement of the circles of radius R centered at the footholds, 
X]"o=i ^"^0 ~ "^^^ above considerations yield the following theorem: 

Theorem 4.10. We can compute 9(JF) n A and the labels of the edges of d{T) 
incident to the arcs of d{J^) Ci A in 0{\A\ logn) time using 0{\A\a{n)) space. 



4.4. Computation of the arcs of d{!F) issued from a foothold. The pre- 
vious section has shown how to compute all the vertices of J- that are incident to at 
least one circular arc. It remains to find the vertices of T incident to two straight 
edges. As we have seen in Section |^, a vertex of T incident to two straight edges of 
d{T) is a foothold. Furthermore, considering a foothold Sig in a cell F of A, Sig is a 
vertex of T incident to two straight edges of d{!F) if and only if Sig is a vertex of the 
convex hull of the footholds reachable from Sig. The k[^ footholds contained in the 
disk of radius R centered at Sig can be found in 0{k[^^) amortized time because we 
have already computed the Delaunay triangulation of the footholds |l^. Thus, we 
can decide if Sig is a vertex of the convex hull of these k[^ footholds in 0{k[^) time 
and space. When Sig is a vertex of the convex hull, we can also find the two edges 
of the convex hull adjacent to Sig in 0{k[^) time and space. As the sum of the k[ for 
i € {1, . . . , n} is bounded by the size of A, we obtain the following theorem: 
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Theorem 4.11. The footholds belonging to d{!F) and the labels of the arcs of 
d{!F) issued from these footholds can be found in 0{\A\) time and space. 

4.5. Construction of T. 

Theorem 4.12. The free space of the spider robot can he computed in 0{\A\ logn) 
time using 0{\A\a{n)) space. 



Proof. By Theorem 4.10, we have computed ah the circular arcs of d{J-) and the 
labels of the edges of d{T) incident to them. By Theorem 4.11, we have computed 



all the vertices of d{J-) that are incident to two straight edges of d{!F) and the label 
of these two edges. It remains to sort the vertices of d{J-) that appear on the line 
segments [st, Sj\. We only consider the line segments [s;, Sj\ such that the correspond- 
ing label {i,j) appears during previous computations. Then, we sort the vertices of 
d{!F) that belong to each such relevant line. Since — Q(\A\) (l), sorting all 

these vertices can be done in 0(|yl| logn) time. A complete description of d{J-) then 
follows easily. □ 

5. Generalization to polygonal foothold regions. 

5.1. Introduction and preliminaries. We consider now the case where the 
set of footholds is no longer a set of points but a set S of pairwise disjoint polygonal 
regions bounded by n line segments ei, . . . , e„. Clearly, 5 is a subset of the free space 
T of the spider robot. Let JFg denote the free space of the spider robot using as 
foothold regions only the edges ei, . . . ,e„. Suppose that the spider robot admits a 
stable placement outside S with its feet inside some polygonal footholds; then the 
placement remains stable if it retracts its legs on the boundary of these polygonal 
regions. Hence, T = J-eU S . We show how to compute JFe- 

the results of Section || remains true if the foothold 



As observed in Remark 3 



regions are line segments provided that Ti^ is replaced by Tig- the generalized helicoidal 
volume defined by (see Figure 5.1): 



Tie 



{{P, 0) £IR^ X \ Pe HD{s, 0), s e ej. 



The helicoidal volume associated to a point site Si will be, henceforth, denoted by 




Fig. 5.1. Section ofTie^ by the "plane" Tig. 



Similarly, we define the generalized circle Ce^ as the set of points at distance R 
from Ci. Let Ae denote the arrangement of the n generalized circles Cej, . . . , Ce^. 
Notice that \Ae\ ^ e(n2). 
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Arc of'conchoid 



Arc of ellipse 
Circular arc ■, 



Circular arc 



Arc of ellipse 



Arg,-6,f conchoid 



Circular arc : 



Fig. 5.2. Example of free space Te, for polygonal foothold regions. The polygonal foothold 
regions are shown in dark grey. The other parts of Te, are in light grey. The Ce^ and some arcs of 
conchoid are dashed. All the line segments tou ching the polygons in two points are of length 2R and 
represent the ladder introduced in Section p. q. 



Each arc of the boundary d{Te) of Te is either an arc of Ce^ corresponding to 
a maximal extension of one leg, or an arc corresponding to placements at the limit 
of stability of the spider robot. Similarly to what we did in Section ^, we compute 
first the contribution of each Ce- to d{Te) (Sections ^.2[) . Thereafter, we compute the 
arcs of d{Te) that correspond to placements where the spider robot is at the lim it of 
stability (Section |]^). Finally, wc show how to construct Te (and T) in Section 5.4 . 

Figure [sT^ shows an example of free space Te for polygonal foothold regions. 

5.2. Computation of d{J-e) n Ae- We compute the contribution to d{!Fe) of 
each generalized circle Ce^ in turn. We consider the contribution of Ce^ to d{Te) 
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for some iq e {1, . . . ,n}. Cg^^ is composed of two half-circles and two straight line 
segments. In order to compute the contribution of Cg^^ to d{J-e), we evaluate first the 
contribution of the half-circles and then the contribution of the straight line segments. 
For convenience, we will not compute the contribution of the half-circles to d{!Fe) but 
the contribution of the whole circles. Similarly, we will compute the contribution of 
the whole straight lines supporting the line segments of Cei^ ■ 

Let and s'^^^ denote the two endpoints of the line segment 6^,, , and let and 
Cs[ denote the unit circles centered at Si„ and s[^^ respectively. Let 1^^ and denote 
the two straight line segments of , and Li^ and L'^^ their supporting lines. We 
show how to compute the contributions of Cs^^ and Li^ to d{J-e)', the contributions 
of Cs' and L'^^ can be computed likewise. 

Let Csig — Csig X and Ci^ = Lj x S^. Basically, we compute d{Te) n Cs-^^ and 



Lig , as explained in Section |4.1| , by computing Ui{Tiei HC^. ), Ui^i,, (Tie; nCs^^^ ] 
Ui {Hei n ) and Ui^ig {Ha n Cig ) . The properties of the new regions . —TLei^ ^si^ 
and = H Lig are different though similar to the properties of Zg- — Tisi ^ ^s; 



described in Section 4.2. The analysis of Z^.^ and J^e^ are subdivided into two parts: 



first, we consider the line Di supporting and we examine the regions Zjj^ = Ti^i; fl 
and 3^£); = Til), n Cig where Ti-Di the generalized helicoidal volume induced by 

A: 

Hd, = {{P. e)eM^ xs^ \ Pe hd{s, e), s e A}. 

Then we deduce Zg. (resp. from Zn., Zs- and Zgi, (resp. yoi, J^s; — H Cig 
and ys'.) where Si and are the two endpoints of e^. Thereafter, we compute the 



contribution of Ce^g to d{!Fe) in a way similar to what we did in Section 4.2. The 
following theorem sums up these results: 

Theorem 5.1. We can compute d{J-e) H Ae and the labels of the edges of d{J^e) 
incident to the arcs of d{Te) Ae in 0(|^e|a7(n) logn) time using 0{\Ae\as{n)) 
space. The proof of this theorem, omitted here, is a direct generalization of the 



proof of Theorem 4.1C . Details are given in M or ||1^ 



5.3. Arcs of d{Te) corresponding to the placements where the spider 
robot is at the limit of stabiUty. We now have to compute the edges of J-g that 
do not belong to Ae- The arcs of d{Te) H Ae correspond to placements at the limit of 
accessibility of the spider robot, and vice versa. Thus, other edges of J^e correspond to 
placements at the limit of stability of the spider robot. We denote by d{!Fe)stab the set 
of those edges. A placement P of the spider robot is at the limit of stability if and only 
if there exists a closed half-disk of radius R centered at P that does not contain any 
foothold except at least two footholds located on the diameter of the half-disk such 
that P is between these footholds (see Figure 5.3). Therefore, the edges of d{J-e)stab 



are portions of the curves drawn by the midpoint of a ladder of length 2R moving by 
translation and rotation such that the ladder touches the boundary of the foothold 
regions in two points but does not intersect the interior of the foothold regions. Hence, 
the edges of d{Te)stab are supported by the projection (onto IB?) of the edges of the 
boundary of the free space of the ladder moving by translation and rotation amidst 
the foothold regions considered as obstacles, i.e., the set of (P, 9) E x IR/ttZ such 
that the ladder of length 2R that has its midpoint at P and makes an angle with 
the X-axis does not properly intersect the interior of the foothold regions. According 
to [0 , the edges of the boundary of the free space of the ladder can be computed in 
0{\Ae\\ogn) time using 0(|^e|) space. The projection (onto IR^) of each edge can 
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easily be computed in constant time. Thus, we can compute, in 0(|^e|logn) time 
and 0(|yle|) space (using ||l6|), a set of curves in Ft^ that contains the arcs of d{!Fe) 
that correspond to placements at the limit of stability of the spider robot. However, 
it remains to compute the portions of these curves that belong to d{J-e). 



Fig. 5.3. Example of placement P at the limit of stability. 



5.3.1. Notations and definitions. The relative interior of an d is called a 
wall. An endpoint of an is called a corner (when several walls share an endpoint, 
we define only one corner at that point) . The ladder is a line segment of length 2R. A 
placement of the ladder is a pair (P, 9) G x M/ttZ where P is the location of the 
midpoint of the ladder and 9 is the angle between the x-axis and the ladder. A free 
placement of the ladder is a placement where the ladder does not properly intersect 
the walls or partially lies on some walls and does not properly intersect the others 
(if none of the polygonal regions of S are degenerated into line segments or points, 
then a free placement of the ladder is a placement where the ladder does not intersect 
the interior of the polygonal regions of 5). A placement of type corner-ladder is a 
placement of the ladder such that the relative interior of the ladder touches a corner. 
A placement of type wall- endpoint is a placement of the ladder such that an endpoint 
of the ladder touches a wall. A placement of type corner- endpoint is a placement of the 
ladder such that an endpoint of the ladder touches a corner. We now define k-contact 
placements of the ladder. 




(a) (b) (c) (d) 

Fig. 5.4. Examples of 2- contact placements of type (a): (corner-ladder)^, (b): (corner-ladder, 
wall- endpoint), (c): (wall-endpoint)^ and (d): (comer- endpoint). 



A 1-contact placement is a free placement of type corner-ladder or wall-endpoint. 
A 2- contact placement is either the combination of two 1-contact placements or a 
free placement of type corner-endpoint. A 2-contact placement is said to be of 
type (corner-ladder)^, (corner-ladder, wall-endpoint), (wall-endpoint)^, or (corner- 
endpoint), in acc orda nce to the types of placements involved in the 2-contact place- 
ment (see Figure 5.4). Given two walls (resp. a wall and a corner, two corners, one 
corner), the set of 2-contact placements induced by these two walls (resp. the wall 
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and the corner, the two corners, the single corner) is caUed a 2-contact curve. The 
type of a 2-contact curve is the type of the 2-contact placement defining the curve. 
Note that the 2-contact curves are defined in IR^ x IR/ttZ. A 3-contact placement 
is a combination of a 1-contact placement and a 2-contact placement. The types of 
3-contact placements are naturally given by (corner-ladder)'^, (corner-endpoint, wall- 
endpoint). . . With this definition, we unfortunately cannot guarantee that all the 2- 
contact curves end at 3-contact placements. Indeed, a 2-contact curve defined by the 
ladder sliding along a wall (see Figure ^^b) ends on one side (if no other wall blocks 
the sliding motion) at a 2-contact placement of type (corner-endpoint), where the 
ladder is coUinear with the wall, without properly intersecting it. In order to ensure 
that all the 2-contact curves end at 3-contact placements, we consider these 2-contact 
placements as 3-contact placements, and denote their type by (corner-endpoint, ||). A 
k-contact placement, > 3, is the combination of p 1-contact placements, q 2-contact 
placements and r 3-contact placements such that p + 2q + 3r — k. 

Now, we define a 2-contact tracing as the projection onto IR^ of a 2-contact curve. 
Similarly as above, we define the types of the 2-contact tracings. Notice that, to any 
point P on a given 2-contact tracing /C, corresponds a unique placement (P, 9) on the 
2-contact curve that projects onto /C. It follows that, to any point P on a 2-contact 
tracing /C, corresponds a unique pair (M, N) of points of contact between the ladder 
at (P, 6) and the walls {M and N are equal when /C is a 2-contact tracing of type 
(corner-endpoint)); when P is an endpoint of /C, a 3-contact placement corresponds 
to P, however, (M, A'') is uniquely defined by continuity. The points M and A^ are 
called the contact points corresponding to P G /C. We also define the three contact 
points corresponding to a 3-contact placement. 

A 2-contact tracing is either a straight line segment, an arc of ellipse, an arc of 



conchoid or a circular arc. Indeed (see Figures 5^, 5_^, 5/7 and 5^), a 2-contact tracing 
of type (corner-endpoint) is a circular arc; a 2-contact tracing of type (wall-endpoint)^ 
is an arc of ellipse; a 2-contact tracing of type (corner-ladder, wall-endpoint) is an 
arc of conchoid (see Q]); a 2-contact tracing of type (corner-ladder)^ is a straight 
line segment. As we said before, we can compute all these 2-contact tracings in 



0(|^e| logn) time using 0(|y^e|) space |16 , and it remains to compute the portions 
of these curves that belong to d{Te). 

5.3.2. Overvievir. We first show that only some portions of the 2-contact trac- 
ings correspond to positions at the hmit of stability of the spider robot (Section ^.3.3| ). 
These portions are called the relevant 2-contact tracings. Then, we prove that we do 
not have to take into consideration the intersections between the relative interior of 
relevant 2-contact tracings (Proposition 5.2). We also show that, if a point A is an 
endpoint of several relevant 2-contact tracings, only two of them can support edges 
of d{Te)stab in the neighborhood of A (Propositions |5.3D . Finally (Section [5.3.4 ), we 
compute a graph whose edges are relevant 2-contact tracings and where the degree of 
each node is at most two. This graph induces a set A of curves supporting d{Te)stab 
(Theorem 5.4) that will allow us to compute d{J-e)stab in Section 5.4. 

5.3.3. Relevant 2-contact tracings. As mentioned above, a placement P of 
the spider robot is at the limit of stability if and only if there exists a closed half-disk 
of radius R centered at P that docs not contain any foothold except at least two 
footholds located on the diameter of the half-disk, one on each side of P. Thus, a 
point P of a 2-contact tracing JC belongs to an arc of d{!Fe)stab only if P lies between 
the two contact points corresponding to P e /C. The portions of the 2-contact tracings 
for which that property holds are called the relevant 2-contact tracings. The other 
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portions are called the irrelevant 2- contact tracings. We now show how to compute 
the relevant 2-contact tracings for each type of contact. Let JC denote a 2-contact 
tracing, let P G J C and le t M an d TV" be the two contact points corresponding to 
P e /C. In Figures 5.5, 5.6, 5.7 and 5.8, the walls and the relevant 2-contact tracings 
are thick, the irrelevant 2-contact tracings are dashed thick, and the ladder is thin. 



Type (corner-endpoint): 

/C is a circular arc, M and N coincide with 
one endpoint of the ladder. Thus, all the 2- 
contact tracings of type (corner-endpoint) are 
wholly irrelevant. 



Type (wall-endpoint)-^: 

K. is an arc of ellipse, M and N are the end- 
points of the ladder and thus, P lies between 
them. Therefore, all the 2-contact tracings of 
type (wall-endpoint)^ are wholly relevant. 



Type (corner- ladder, wall-endpoint): 



JC is an arc of conchoid. If the distance be- 
tween the corner and the wall is greater than 
R, then JC is wholly relevant. 



Otherwise, if that distance is smaller than R, 
then, the two relevant portions and the irrel- 
evant portion of JC are incident to the corner 
involved in the type of JC. 

Notice that, if t he c orner is an endpoint of the 
wall (see Figure 5^b), then IC degenerates into 
a line segment and the irrelevant portion of IC 
is the portion which is not supported by the 
wall. 

Type (corner- ladder)'^: 

/C is a line segment. If the distance between 
the two corners is greater than _R, then JC is 
wholly relevant; otherwise, the portion of JC 
which is relevant, is the line segment joining 
the two corners. 



Fig. 5.5. Irrelevant 2-contact tracing of 
type (corner-endpoint), i.e., circular arc. 




Fig. 5.6. Relevant 2-contact tracing of type 
(wall-endpoint)'^ , i.e., arc of ellipse. 




Fig. 5.7. Relevant, and partially relevant, 
2-contact tracings of type (corner-ladder, 
wall-endpoint), i.e., arcs of conchoid. 




Fig. 5.8. Relevant, and partially relevant, 
2-contact tracings of type (corner-ladder)-^ . 
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We now show that the intersections between the relative interiors of the relevant 
2-contact tracings are not interesting for the spider robot motion problem. We recall 
that, if a vertex A of d{J-e) belongs to Ae, then we know by Theorem 5.1 the labels of 
the edges of d{J-e) incident to A. Otherwise, if A ^ Ae, then the two edges of d{Te) 
that end at A correspond to placements at the limit of stability of the spider robot. 

Proposition 5.2. Any vertex A of d{J^e), such that A ^ Ae, is an endpoint of 
the two relevant 2-contact tracings supporting the edges of d{!Fe) ending at A. 

Proof. Since the two edges of d{J-e) that end at A correspond to placements at 
the limit of stability of the spider robot, they are both supported by some relevant 
2-contact tracings. Thus, we only have to prove that A is an endpoint of these two 
relevant 2-contact tracings. 

Let /Ci and /C2 be these two relevant 2-contact tracings and assume for a con- 
tradiction that A is not an endpoint of /Ci (nothing is assumed for A with respect 
to /C2)- Let Li = {A,9i) (resp. L2 = {A,02)) be the placement of the ladder that 
correspond to ^ e /Ci (resp. A G JC2) and let Mi and A^i (resp. M2 and be the 
corresponding contact points (see Figure 5.9). First, notice that Li ^ L2. Indeed, 
otherwise, Li is at least a 3-contact placement and then, A must be an endpoint of 
/Ci, which contradicts our assumption. 




Fig. 5.9. For the proof of Proposition 5. 



By the definition of the relevant 2-contact tracings, A is between Mi and A^i. 
Moreover, A cannot be equal to Mi or A^i since A is not an endpoint of JCi. It follows 
that neither M2 nor N2 is equal to A, because otherwise Li would be a 3-contact 
placement. Therefore, A is strictly between Mi and A^i, and strictly between M2 and 
^^2- Thus, A is strictly inside the polygon {MiM2NiN2). 

On the other hand, since A ^ Ae, A does not belong to any Ce^, and therefore, 
the walls supporting A/i, A^i, M2 and N2 intersect the open disk Da of radius R 
centered at A. Thus, there exists four points M(, iV(, M'2 and N'2 on these walls and 
in Da, that are close enough to A/i, A^i, M2 and N2 respectively to ensure that A 
belongs to the interior of the polygon (M(M2iV( A^j)- Since the distances from A to 
M{, iV(, M[ and iVj, are strictly smaller than _R, A belongs to the interior of Te- This 
contradicts our assumption that A is a vertex of d{Te) and yields the result. □ 

Consider now the adjacency graph Q of the relevant 2-contact tracings such that 
two relevant 2-contact tracings are connected in Q if and only if they have a common 
endpoint (the intersections between the relative interiors of the relevant 2-contact 
tracings are not considered). Notice that, given the set of relevant 2-contact tracings, 
Q can be easily computed in 0(|yle| logn) time. Now, given two vertices of d{!Fe)C^Ae 
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that are connected along d{J-e) by arcs of d{J-'e)stabi we want to compute these arcs. 
For computing these arcs, we cannot simply use the graph Q because the degree of 
some nodes of G may be arbitrarily large (see Figure 5.10| ). We show in the next 
proposition that we can deduce from Q a graph Q* such that the degree of each 
node of G* is at most two and that G* supports any portion of d{!Fe) which is the 
concatenation of arcs of d{!Fe)stab- 



- ^^^^ 
















Fig. 5.10. Relevant 2-contact tracings K,i, . . . , IC7 ending at A. K\,K.2,K,3 and IC5 are 2- 
contact tracings of type (corner-ladder, wall-endpoint) (i.e., arcs of conchoid). is a degenerated 
2-contact tracing of type (corner-ladder, wall-endpoint) (i.e., a line segment). Ka and K5 are 2- 
contact tracings of type (corner-ladder)^ (i.e., line segments). 

We consider four hypotheses (HI,. . .,H4) that obviate the need to consider degen- 
erate cases. They are not essential but substantially simplify the proof of the following 
proposition. The first three hypotheses are made to ensure that the degree of each 
vertex of the free space of the ladder is three. 

HI The line segments ei, . . . , e„ compose the boundary of a set of non degenerated 

polygons (i.e., no polygon is reduced to a line segment or to a point). 
H2 The ladder does not admit any 4-contact placement. 

H3 The arc (of conchoid) drawn by an endpoint of the ladder when its other endpoint 
moves along a wall while the ladder remains in contact with a corner, is not 
tangent to any other wall. 

H4 The ladder does not admit any 3-contact placement when its midpoint is located 
at a corner. 

Proposition 5.3. For any node A of G of degree k such that A ^ Ae, at most 
two relevant 2-contact tracings can support d{(Fe) in a sufficiently small neighborhood 
of A. Moreover, we can determine these at most two curves in O(fclogfc) time using 
0{k) space. 

Proof. Let A ^ Ae he a, node of G of degree k. We assume that k > 2, otherwise 



Proposition 5.3 is trivial. Let /Ci, . . . ,/Cfc be the relevant 2-contact tracings that end 
at A, and let Li = {A, (j>i) be the placement of the ladder that corresponds to ^ e JCi. 
Da is the open disk of radius R centered at A. We distinguish two cases whether A 
is a corner or not. 



Case 1: ^ is a corner. (See Figure 5.10 



The 2-contact tracing /C^ involves at least another contact than the corner-ladder 
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Fig. 5.11. Wedge Pi-^APi^ is in near A. 



contact at A. This contact cannot be of type corncr-cndpoint by Hypothesis H4. If 
the contact is of type wall-endpoint, we define Pi as the contact point between this 
wah and the ladder at placement Li (see Figure 5.11 ). Since A ^ Ae, the wall must 
intersect Da and we define P/ as a point close to P; in that intersection. If the contact 
is of type corner-ladder, we define P; = P/ as the corner (distinct from A) involved in 
this contact (notice that P; = P/ e Da by Hypothesis H4). 

Fact: Vi ^ j, (j)i ^ 
Otherwise, Li = Lj is a 3-contact placement contradicting Hypothesis H4. 

Fact: A is a non-flat vertex of CH{A, Pi, . . . , Pfc) or belongs to the interior of Te- 
Assume that A g d{!Fe). Then, A lies on the boundary of CH(A, Pi, . . . , Pk) because, 
otherwise, the P/ provide footholds such that the spider robot can move in a neigh- 
borhood of A. Furthermore, A must be a non-flat vertex of CH(A, Pi, . . . ,Pfc), by 
Hypothesis H4. 

Assume now that A £ d{J-e), and let P^^ and P^^ be the two vertices of CH(A, 
Pi, ... , Pk) such that Pi^ , A and Pi^ are consecutive along the boundary of CH(A, 



Pi, ... , Pfc) (sec Figure 5.11 ). We will exhibit a stable placement for the spider robot 
at any position P inside the intersection of the wedge P^ ^ APi,^ and a neighborhood of 
A. Let hi and ft,2 be two points in the wedge Pi-^APi^ such that the wedges Pi^Ahi 
and h2APi^ are right (see Figure 5.11). 

— If P is in the wedge Pi^Ah2, and is close enough to A, the footholds A, Pi^ and 
P/^ yield a stable placement for the spider robot. 

— If P is in the wedge h2Ahi, and is close enough to A, footholds A, P[^ and P/^ 
yield a stable placement for the spider robot. 

— If P is in the wedge hiAPi-^, and is close enough to A, footholds A, P/^ and Pi^ 
yields a stable placement for the spider robot. 

Fact: ICi, i ^ {11,12}, cannot support an edge of d{!Fe) incident to A. 
We assume that A E d{!Fe) because, otherwise, the claim is obvious. It follows that A 
is a non-flat vertex of CH(A, Pi, ... , Pk). A 2-contact tracing /C^, i ^ {ii, 12}, cannot 
be an arc of ellipse because, otherwise, Li is a 3-contact placement (because A is 
here a corner) contradicting Hypothesis H4. Then, JCi can be either the segment APi 
or an arc of conchoid. If K,i is an arc of conchoid, then, by the general properties of 
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conchoids (see [||), /C^ is tangent to the segment APi at A. Thus, JCi is always tangent 
to the segment APi at A. The point Pi strictly belongs to the wedge Pi^APi^ , because 
we have shown that (j)i ^ , c^ij }. Thus, in a neighborhood of A, JCi is strictly inside 
the wedge Pi^^APi^ and thus strictly inside J^^- Therefore, )Ci cannot support d{J-e), 
in a neighborhood of A. 

Hence, by sorting the Pi by their polar angles around A, we can determine, in 
O(fclogfc) time, if A is a non-flat vertex of CH(A, Pi, . . . ,Pk), and if so, determine 
ii and i2. If A is a non-flat vertex of CH(A, Pi, ... , P^), then, only /C^^ and JCi^ can 
support an edge of d{Te) incident to A. Otherwise, A belongs to the interior of 
and none of the 2-contact tracings /Ci, . . . ,/Cfe can support an edge of incident 
to A. 

Case 2: A is not a corner. 

Fact: // there exists i ^ j such that (pi ^ (pj, then A belongs to the interior of 

For each relevant 2-contact placement Li — (A, (/j^), there exists two contact points Mj 
and Ni on each side of A at distance less or equal to R. Since A is not a corner, neither 
Mi nor Ni is equal to A, thus A belongs to the relative interior of the segment MiNi. 
It follows, when (pi ^ pj, that A belongs to the interior of the polygon (MiMjNiNj) 



(see Figure 5.9). Similarly as in the proof of Proposition 5.2, since A ^ Ae, there exists 
four footholds M'i, N'i, M'^, N'^ in Da and in some neighborhoods of Mi,N,Mj,Nj, 
respectively, such that A belongs to the interior of the polygon (MiMjNiNj). Thus, 
A belongs to the interior of J-e- 

Hence, if there exists i ^ j such that (j)i ^ pj, none of the 2-contact tracings 
/Ci, . . . , /Cfc can support an edge of d{!Fe) incident to A. We now assume that pi = (pj, 
Vi,j. 

Fact: There are at most six 2-contact tracings incident to A. 
The general position hypothesis H2 forbid fc-contacts for fc > 3, thus A corresponds 
to a 3-contact placement. The three possible choices of two contacts among three, 
give three 2-contact tracing intersecting in A and thus, six arcs incident to A. 

Fact: There are three 2-contact tracings incident to A. 
If the 3-contact placement L is of type (corner-endpoint, ||), then there are only three 

2- contact tracings incident to A, that are two circular arcs and one line segment. 
Otherwise, it comes from the general position hypotheses HI, H2 and H3 (designed 
to ensure that property) that a 2-contact tracing cannot be valid on both side of the 

3- contact, i.e., on one side of the 3-contact placement, the placements are not free. 
The proof that the hypotheses ensured that fact is detailed in Q . 

Fact: There are two relevant 2-contact tracings incident to A. 
Since A is not a corner, at the 3-contact placement i, two contact points are on the 
same side of A. Thus, only two of the three 2-contact tracings incident to A are 
relevant. □ 

5.3.4. Construction of A. Now, consider the graph Q and each node A in turn. 
If j4 e yle, we disconnect all the edges of Q that end at A. Notice that for each such 
node A^ we know, by Theorem whether A € d{Te) and, in such a case, the labels 
of the edges of d{Te) incident to A. li A ^ Ae, we disconnect the edges ending at 
A except th ose (at most two) that may support d{Te) in a neighborhood of A (see 



Proposition |5.3| ). In this way, we obtain a graph Q* such that the degree of each node 
is one or two. We consider each connected component of this new graph as a curve. 
Let A be this set of curves. These curves are represented in Q* as chains (open or 
closed). It follows that, even if a curve is not simple, there exists a natural order 
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along the curve. Then, according to Propositions 5.2 and |5.3| , we get the following 
theorem: 

Theorem 5.4. We can compute, in 0{\Ae\ logn) time using 0{\Ae\) space, a set 
A of curves that support the edges of d{!Fe) corresponding to placements at the limit 
of stability of the spider robot. Moreover, any portion V of d{!Fe) either intersects Ae 
or belongs to a unique curve of IS.. 

5.4. Construction of Te and T. We can now construct Te and T. Let Xk{n) 
denote the maximum length of the Davenport-Schinzel sequence of order A; on n 
symbols and ak{n) = Xk{n)/n. Note that a^ln) — a{n). 

Theorem 5.5. Given, as foothold regions, a set of n non intersecting straight 
line segments that satisfies Hypotheses HI, H2, H3 and H4-, we can compute the free 
space Te of the spider robot in 0(\Ae\cii,{n)^0'gri) time using OQAelasin)) space. 

Proof. By Theorem |5.l| , we can compute the contribution of Ae to d{Te) and 
the label of the edges of d{Te) incident to them in 0(|.4e|a7(n) logn) time using 
0(|^e|a8('^)) space. By Theorem |5.4 we can compute, in 0{\Ae\\ogn) time using 
0(|vAe|) space, a set A of curves that support the edges of d{J-e) that do not belong 
to Ae- Moreover, any portion V of d{!Fe) such that V H Ae =0 belongs to a unique 
curve of A. Thus, by sorting all the vertices of d{Te) n .4e n A on the relevant curves 
of A, we obtain all the edges of d{!Fe) that belong to a connected component of 
d{!Fe) intersecting Ae- Indeed, for each vertex A S d{J-e) n .4e n A, we know, in 
a neighborhood of A, the portion of the curve of A that belongs to d{Te) because 
we can simply determine, for each edge, a side of the edge that belongs to Te (the 
contact points corresponding to the edges determine a side that necessarily belongs 
to Te^. Then, it is an easy task to deduce all the connected components of d{!Fe) 
that intersect Ae- 

It remains to compute the connected components of d{!Fe) that do not intersect 
Ae. Each of these components must be a closed curve of A. Moreover, all the curves 
of A belong to Te. Thus, according to Theorem |S.4| , any closed curve /C of A that 
does not intersect Ae is either a connected component of d{Te) or is strictly included 
in Te. Therefore, by considering, in addition, all the closed curves of A that do not 
intersect Ae, we finally obtain a set 5* of closed curves that contains d{Te) and such 
that any curve of ^' is either a connected component of d{Te) or is strictly included 
in Te. 

At last, as we can simply determine, for each curve of a side of the edge 
that belongs to Te, we can easily deduce from ^ the free space Te. That concludes 
the proof since all these computations can be done in 0(|.4e|a8('T') log ft) time using 
0{\Ae\as{n)) space. □ 

As we said at the beginning of Section |^, the free space of the spider robot using 
as foothold regions a set of polygonal regions is obtained by adding these polygonal 
regions to Te- This does not increase the geometric complexity of the free space nor 
the complexity of the computation. Thus, we get the following theorem: 

Theorem 5.6. Given a set of pairwise disjoint polygonal foothold regions with n 
edges in total that satisfies Hypotheses HI, H2, H3 and H4, we can compute the free 
space T of the spider robot in 0{\Ae\cta{n) logn) time using 0{\Ae\cxs{n)) space. 

The function a8{n) is extremely slowly growing and can be considered as a small 
constant in practical situations. This result is almost optimal since, as shown in 
r2(|.4e|) is a lower bound for the size of T. 



^Observe that when the edge belongs to J^e, its two sides belong to J^e- 
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6. Conclusion. We have seen in Theorem 4.12 that, when the foothold regions 
are n points in the plane, the free space of the spider robot can be computed in 
0(|^|logn) time using 0(|yl|a(n)) space where a{n) is the pseudo inverse of the 
Ackerman's function and A the arrangement of the n circles of radius R centered at 
the footholds. By § the size of J- is known to be G(|.A|). The size of A is O(n^) but 
it has been shown in that |^| = 0{kn), where k denotes the maximum number of 
disks of radius R centered at the footholds that can cover a point of the plane. Thus, 
in case of sparse footholds, the sizes of A and J- arc linearly related to the number of 
footholds. 

When the foothold regions are polygons with n edges in total, the free space of the 
spider robot can be computed in 0(|^e|a8(n) logn) time using 0{\Ae\as{n.)) space, 
where nak{n) = Afc(n) is the maximum length of a Davenport-Schinzel sequence of 
order fc on n symbols, and Ae is the arrangement of the n curves consisting of the 
points lying at distance R from the straight line edges. Note that the size of Ae is 
0(n2). 

It should be observed that, in the case of point footholds, our algorithm im- 
plies that 0(|^|a(n)) is an upper bound for \J-\. However, this bound is not tight 
since |JF| = 0(|y^|) [y. In the case of polygonal footholds, our analysis implies that 
0(|yie|Q!8(?T-)) is an upper bound for \!F\. We leave as an open problem to close the 
(small) gap between this upper bound and the r2(|^e|) lower bound. 

Once the free space is known, several questions can be answered. In particular, 
given two points in the same connected component of the algorithm in [Q computes 
a motion of the spider robot, i.e., a motion of the body and a corresponding sequence 
of leg assignments that allows the robot to move from one point to the other. 

The motion planning problem for other types of legged robots remains to be 
studied. The case where all the legs are not attached at the same point on a polygo- 
nal/polyhedral body is particularly relevant. A spider robot for which all the legs are 
not of the same length is also an interesting model. 

Acknowledgments. We would like to thank Joseph O'Rourke for helpful com- 
ments. 
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